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ABSTRACT 

The high resolution image of the jet of the BL Lac object MKN501 in radio, show a 
limb-brightened feature. An explanation of this feature as an outcome of differential 
Doppler boosting of jet spine and jet boundary due to transverse velocity structure 
of the jet requires large viewing angle. However this inference contradicts with the 
constraints derived from the high energy 7-ray studies unless the jets bends over a 
large angle immediately after the 7-ray zone (close to the central engine). In this 
letter we propose an alternate explanation to the limb-brightened feature of MKN501 
by considering the diffusion of electrons accelerated at the boundary shear layer into 
the jet medium and this consideration does not require large viewing angle. Also the 
observed difference in the spectral index at the jet boundary and jet spine can be 
understood within the frame work of shear acceleration. 

Key words: galaxies: active - galaxies: jets - BL Lacertae objects: individ- 
ual(MKN501) - acceleration of particles - diffusion 



1 INTRODUCTION 

BL Lac objects are the extreme class of active galactic nu- 
clei(AGN) with weak or no emission lines and are cate- 
gorized along with flat spectrum radio quasars(FSRQ) as 
blazars. Their spectra cover a broad range of photon ener- 
gies starting from radio to gamma rays with a few of them 
detected in TeV energies by gr o und based Air Cerenkov 



experiments dKrawczvnskil (12004 ): Katarzvhski. Sol, fc Kusl 
l|200ll) : ISambrunal (|2000l ) lcostamante fc Ghisellinil l|2002)). 



These sources are found to be strongly variable with 
flare time scales ra nging from days to less than an 



hour dGaidos et all dl996l); ICoppi fc Aharonian] l|l999l ) 



ISambrunal i|2000l ): iKrawczvnski et all l|2000l) ). The short 



time variability and their detection at very high energies 
demand that the emission region should be moving down a 
jet at r elativistic velocit i es clos e to the line of sight o f the o b- 
server l|Ghisellini et all (|l993l ): iDondi fc Ghisellinil (119951) 1. 
The strong polarization detected in radio/optical energies 
and the non-thermal photon spectra indicates the radio to 
x-ray spectra is due to synchrotron radiation from a non- 
thermal electron distribution cooling in a magnetic field. 
However the gamma ray emission from these sources is 
still not well understood. Leptonic models explain the high 
emission as inverse Compton scattered synchrotron photons 
by the electron population responsible for the synchrotron 
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process itseff(SSC) dMaraschi, Ghisellini. fc Celottil (|l992l) : 
iBloom fc Marscherl l|l996l ): iBottcherl |2000l )) where as in 
hadronic models it is due to the synchrotron proton emission 
and proton-photon interactions involving an exte rnal photon 
field ( s ynchrotron proton blazar model(SPB)) jMannheiml 
1 19981 ) ; lMucke et al l (|2003l )). Under unification hypothesis of 
radio-loud AGN, BL Lac objects are considered to be aligned 
j et version of Fanaroff- R iley type I (FRI) radio galaxies 
dUrrv fc Padovanil l|l995l )). 

MKN501 is a nearby BL Lac object (z=0.034) and also 
the second extra galactic source detected in Te V photon en- 
ergies by ground based Cherenkov Telescopes (|Quinn et alj 
(1996)). It was later detected in Me V photon energies by t he 
satellite based experiment EGRET (|Kataoka et a l. (1999)). 
The radio images of MKN501 show a j et emerging from 
a brig h t nucleus dEdwards et al" d2000h : Giovannini et al] 
1 19991 ): lAaronl dl999l ): IGiroletti et all d2004l )). The high 
resolution (milli arc second) radio images show a trans- 
verse jet structure with the edges being brighter than 
the centra l spine commonly ref er red as "limb-brig h tened " 
structure dEdwards et al.l d2000l ); iGiovannini et alj (|l999h ; 
IGiroletti et al.l ( 2004 )). This feature is usually explained 
by the "spine-sheath" model where the velocity at the jet 
spine is larger compared to the velocity at the boundary. 
Such a radial stratification of velocity across the jet arises 
when jet moves through the ambient medium and the vis- 
cosity involved will cause a shear at the boundary. Three- 
dimensional hydrodynamic simulations of relativistic jets 
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l|Alov et all (|2000l )) and two-dimensional simulations of rel- 
ativistic magnetized jets l|Leismann et al.l l|2005l )) also sup- 
ports the presence of jet velocity stratification due to its 
interaction with the ambient medium. The existence of ve- 
locity shear at the jet boundar y was first suggested by 
lOwen. Hardee, fc Cornwelll (1 19891) to e xplain the morphol- 
ogy of M87 jet. Perlman et alj (fl999'l later confirmed it 
through the polarisation studies of M87 jet. If the jet is mis- 
aligned towards the observer, it may happen for a proper 
combination of velocities we see a Doppler boosted image 
of the boundary compared to the less boosted sp i ne giv - 
ing rise to a limb-brightened structure (|Komissaro 
iLaind (|l996l )). A possible consequence of the velocity shear 
is the alignment of the magnetic field at the boundary par- 
allel to the flow velocit y due t o stre tching of the frozen-in 
field lines of the plasma l|Kahnl |l983l) ). The polarisation an- 
gle observed at the jet boundary of MKN5 1 is perpendi c- 
ular to the jet axis l|Pushkarev et all i|2005l ); lAaronl (|l999h ) 
indicating a parallel magnetic field. However it should be 
noted here that the polarisation angle at the jet spine in- 
dicates a perpendicular magnetic field and this along with 
the parallel magnetic field at the jet boundary can be 
an outcome of a dynamically dominan t toroidal magneti c 
field structure (iPushkarev et all (l2005h : iGabuzdal [| 19991) : 
iGabuzda, Murray, fc Croninl ( 20051 )). The radial velocity 
stratification of the jet can introduce Kelvin-Helmholtz in- 
stability and the stability of j ets against this in st abilit y 
was studied by many aut h ors ITTurland fc Scheuei (I 19761); 
Blandford fc PringleH 19761); iFerrari. Trussoni. fc Zaninett: 
il 19781) jlHardeel (ll97d);lBirkmshawl i|l99ll )). 

iGiroletti et al.l (|2004l ) have studied the limb-brightened 
structure of MKN501 jet considering the differ ential Dopple r 
boosting at t he je t spine and the boundary (|Laind (1996); 
iKomissarovl l|l990l )) and concluded the viewing angle (an- 
gle between the jet and the line of sight of the ob- 
server) of the radio jet should be more than 15°. How- 
ever high energy studies of MKN501 demands the view- 
ing angle of the jet should be m 5° in order to ex- 
plain th e observed rapid variabil it y and the high energy 
emissi on (|Katarzvriski. Sol, fc Kusl |200ll) ; iTavecchio et al.l 
(2001)). Considering the fact that the gamma ray emission 
is originated fr o m the inner part of the jet close to nucleus, 
IGiroletti et al.l (|2004l ) suggested a bending of the jet may 
happen immediately after the gamma-ray zone to explain 
the required large viewing angle of the radio jet. However the 
mechanism required to bend the jets are still not well under- 
stood (jets deflected due to the pressure gra d ient in external 
medium is studied bvlCanto fc Ragal |l996); Rag a fc Cantol 
l|l996h ; iMendoza fc Longairl (|200ll )) and moreover the ob- 
served large bending of the jet in the radio maps can be 
apparent one because of projection effects. This projection 
effects are even amplified when the jet is close to the line of 
sight. Though it needs to be note d here that jets w i th lar ge 
bending angle are indeed observed (|Savolainen et all |2006)). 

The limb-brightened structure can also be explained 
if we consider the synchrotron emission from the par- 
ticles accelerated at the boundary and this in f erence 
does not require large viewing angle. lEilekl (|l979l . I1982T ) 
considered the acceleration of particles due to tur- 
bulence initiated by Kelvin-Helmholtz and Rayleigh- 
Taylor instabilities at the jet boundary. Particles at the 
boundary can also be accelerated via shear accelera- 



tion ( Bcrezhko (|l98ll ); iBerezhko fc Krvmskiil i|l98ll )) and 
this case is considered in the present work. The accel- 
eration of particles in a shear flow or by turbulence 
is well studied by various authors for both rel a tivisti c 
and non relativistic case (lEarl. Jokipii. fc Morfil] dl98Sl); 
IWeb ' 
(200 



Stawarz fc Ostrowskil 



.19891); lOstrowskil (ll99cK7 . 

_ Rieger fc Duffvl (l2006l ); IStawarz fc Petrosianl (|2008h ; 

IVirtanen fc Vainiol l|2005l )). 



In this letter we explain the observed limb-brightened 
feature of MKN501 by considering the diffusion of electrons 
accelerated at the jet boundary via shear acceleration. In the 
next section we show the required condition for the shear ac- 
celeration to be dominant over turbulent acceleration and in 
§3 we consider the diffusion of particles accelerated at the 
boundary into the jet medium. In §4 we discuss the spec- 
tral index of the particle distribution accelerated via shear 
acceleration process and turbulent acceleration process and 
show the observed index at the boundary of MKN501 jet 
supports the earlier case. Throughout this work, H = 75 
km s" 1 Mpc -1 and qo = 0.5 are adopted. 



2 SHEAR ACCELERATION AT MKN501 JET 
BOUNDARY 

The particle acceleration process at the jet boundary can be 
described by the diffusion equation in momentum space. The 
ev olution o f an is otropic phase space distribution is given 
bv j|Melrosel (|l968D ) 



df{p) 
dt 



l_ d_ 

p 2 dp 



* nf ,df( P )\ 
PD{P) ^P-) 



(1) 



where D(p) is the momentum diffusion coefficient. The char- 
acteristic acceleration timescale can be written as 



tac 



d_ 

dp 



(p 2 D(p)) 



(2) 



If we consider a sheared flow, the electrons are scattered 
across different velocity layers b y turbulent struc tures which 
are embedded in the shear flow. iBerezhkol (|l98ll ) showed in 
such case there will be a net gain of energy in the elec- 
trons getting scattered and this process is referred as shear 
acceleration. The momentum diff usion coefficient i n cas e 
of a shear flow can be written as feieger fc Duffvl l|2006l ); 
iRieger. Bosch-Ramon, fc Duffvl (|2007r i) 



D s {p) = XV t 



(3) 



where r is the mean scattering time given by r ~ A/c with 
A the mean free pat h and y is the sh e ar coe fficient given for 
a relativistic flow as l|Rieger fc Duffvl l|2004 )) 



X : 



15(r(r) 2 



1 



(4) 



where T(r) is the bulk Lorentz factor of the flow and r is 
the radial coordinate of the jet cross section. Using @, the 
shear acceleration timescale(iicc) for r = r a p^ will be 

f 



(4 + £)xt 



(•») 



In case of turbulent acceleration(stochastic), the particles 
are scattered off by randomly moving scattering centres and 
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gets energized by second order Fermi acceleration. The mo- 
mentum diffusion coefficient in this case can be approxi- 
mated as (|Rieeer. Bosch-Ramon, fc Duffvl (|2007n ) 



Dt(p) 



p_ 

3t 



where the Alfven velocity (Va) is given by 

B 



Va = 



(6) 



(T) 



here B is the magnetic field and p the mass density of the 
jet. Hence the turbulent acceleration timescale(ti t 2 c ) will be 



3t 



(4-0 



(8) 



For shear acceleration to be dominant over turbulent acceler- 



ation ticc < tacc- If we consider Bohm diffusion (£ = 1) then 
the mean free path of the electron align ed to the magnetic 
field (Ay) scales as the gyro radius fr„) (|Achterberg fc Ball 

l|l994l )). A || ~ r\ "'™g C , where r/ is a numerical factor (r\ > 1 
for magnetized particles) and 7(3> 1) is the Lorentz fac- 
tor of the electron scattered. Since the magnetic field at 
the jet b oundar y of M K N501 is parallel t o the je t axis (or 
toroid al) (|Aaronl (|l999l l; iPushkarev et all i|2005l ); IGabuzdal 
(1 19991 )). we consider r ~ Ay/c. Also if we consider 



or 

dr 



AT 

Ar 



(9) 



where AF is the difference between the bulk Lorentz factor 
at the jet spine and the jet boundary and Ar is the thickness 
of the shear layer, then the condition for shear acceleration 
to be dominant over turbulent acceleration will be 



Ar < 



yym e c 3 (AT) 
eB 2 



4irp 



3(r(r)" - 1) 



(10) 



If we consider the mass density of the jet is dominated by 
cold protons and if the number of protons are equal to the 
number of non-thermal electrons, then the jet mass density 
can be written in terms of equipartition magnetic &eld(B eq ) 
as 



m p B 2 q (2a 



1) 



167rm e c 2 a7„ 



and (fTOl) will be 



Ar < 0.29 



777c 2 (Ar) 

eB™ 



e rn p (2a — 1) 



tt7„„(r(r) 2 - 1) 



(11) 



(12) 



is the 



where a is the observed photon spectral index, 
proton mass and jmin is the Lorentz factor of electron re- 
sponsible for the minimum observed photon frequency v m i n . 
The equipartition magnetic field can be expressed in terms 
of observed quantities as 



B e 



9.62——(m e cev 71 
T(r) 



d 2 L F(v„ 



Va T (2a 



G (13) 



where F(v m i n ) is the flux at the minimum observed fre- 
quency Vmin, dh is the luminosity distance, V is the vol- 
ume of the emission region and ot is Thomson cross sec- 
tion. Hence, for T(r) 2 ^ 1 and a ~ 0.7, shear acceleration 



will dominate the particle spectrum at the jet boundary of 
MKN501 if the thickness of the shear layer 



Ar < 7.22 x 10" 



F(10MHz) 



/ Ar 



Vobs 



(- 



1.6GHz/ VlOMHz 



910mJy 



R 



1.5parsec 



parsec 



(14) 



Where R is the radius of the spherical region considered. 
(We assume 10MHz as minimum observed frequency and 
the flux at 10MHz is obtained from the flux at 1.6GHz 
considering the same spectral index. The flux at 1.6GHz and 
R in (|14[) are obtained from a region around R.A 10 mas an d 
declination —10 mas from Fig. 7 of lGiroletti et al.l l|2004) L 
The corresponding equipartition magnetic field B eq for T = 
5 is 1.2 x 10 _3 G. 

The electrons accelerated by shear acceleration cool via 
synchrotron radiation. The cooling time for synchrotron loss 
is given by 

67rm e c 



tcool - 

Using ([5| and dT5j) , we find 
~ 1.5 x 10 

tcool 

Ar 



(15) 



10~ 9 parsec 



1.2 x 10 
J / AT 

To" 



3Q 
2 



(16) 



(s) 

and since t a cc -C t coo i, shear acceleration dominates over 
synchrotron cooling. It can be noted that (|16|) is independent 
of the electron energy and hence the maximum energy of 
the electron will be decided by the loss processes other than 
synchrotron loss (which are not considered in this simplistic 
treatment). 

If we maintain the general form of mean scattering time 
t = rop^ , then for shear acceleration to dominate over turbu- 
lent acceleration the thickness of the shear layer (Ar) should 
be 



Ar < 1.7 x 10 



Top s (Ar) 

T(r) 



(4 + C)(2a- 



a(4 - £)7n 



(17) 



It can be noted that (|10[) is equal to (1 1 71) if we set in the 
latter £ = 1 and rop 5 = r/r 9 /c. 



3 PARTICLE DIFFUSION AT THE JET 

BOUNDARY AND LIMB-BRIGHTENING 

Particles accelerated at the shear layer of the jet bound- 
ary, diffuse into the jet medium before getting cooled off 
via synchrotron radiation. As the magnetic field at the jet 
bound a ry is parallel to the j et axis (or toroid al) (I Aaronl 
l|l999l ): IPushkarev et all l|2005h : IGabuzdal (|l999h ). the ra- 
dial diffusion of the electron into the jet medium is deter- 
mined by cross field diffusion. The cro s s field diffusion coeffi - 
cient can be approx i mate d as l|Axfordl (1 19651 ) ; [jokipiil l|l987l ); 
lAchterberg fc Balll (|l994p ) 



377 



(18) 



Where r/(> 1) is the scaling factor determining the field 
aligned mean free path (see (31])). 
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The radial distance Rdiff that the electron diffuse be- 
fore getting cooled can then be approximated as 



Rdiff ~ V K±t coo l 



(19) 



Using (|15l) and (|18[) and considering the equipartition mag- 
netic field we get 



i?^ ^ 2.9 x lrr 4 (^p (j 



B 



.2 x 10" 3 G* 



parsec 



(20) 

Since the thickness of the shear layer Ar <C Rdi / / (refer 
(|14[1 and (|20p ). the thickness of the limb brightened structure 
will be ~ Rdiff - This corresponds to an angular distance of 
4.7 x W~ 4 rnas which is beyond the resolution of present day 
telescopes. 

For r = Top^ , the cross field diffusion coefficient will be 



A 2 -£ 

3r 



(21) 



Using JT5JI and we get 

J? dl// ~ 5.2 x W^B^To^p^cm (22) 

and hence the thickness of the limb brightened structure will 
be energy dependent for £ ^ 1. 



4 SPECTRAL INDEX 

If we add mono-energetic particle injection term (5(p — p )) 
and particle escape term ( — l/t esc ) in fl]), then the steady 
state equation in case of shear acceleration for p > p and 
£ = 1 can be written as 



dp^ dp 



: 



and in case of turbulent acceleration it will be 

ft 



d 2 ft | dU 
dp 2 dp 



1pt e 







(23) 



(24) 



where ip - 



Zc £ t . If we substitute p = 1/a; in (|23[) we get 



3^ 
dx 







(25) 



dx 2 da; x T <>te 

Equat ions (|24[) and (125[) can be solved analytically l)Kepinskil 
() 1905T ) ^ and the solutions are complex and are given by 

/.= ' 1 

X 

and 

ft = 



X T oPtesc 

a s J4 ( 2i 



Tpt e 



XToPU 



atJi 2i 



+ b s Yi 2i 



1 



X T oPte 



(26) 



+ b t Y 2 2i 



(27) 

Where J n {z) and Y n (z) are the Bessel functions of first and 
second kind and a s , b s , a t and bt are constants. For negli- 
gible esca pe (t esc — > oo), using the limi ting forms of Bessel 
functions (|Abramowitz fe Stegunl l|l972t )1. the solutions (|26[) 
and (|27p approaches a power law / s oc p~ 4 and ft oc p~ 2 . 
The shear accelerated particle number density will then be 
n s (p) oc p~ 2 and the corresponding synchrotron photon flux 



will be Sv l3 hear oc v~^^ 2 . For turbulent acceleration the 
number density will be independent of p («t(f>) oc p°) and 
hence the observed synchrotron photon flux will be a flat 
one Sti turb 

(xv'. The spectral index map of MKN501 jet 
indicates a steep ph oton spectra at the bo undary and a flat 
spectra at the spine (|Giroletti et al.l l|2004h ). Hence it can be 
argued that the shear acceleration may be dominant at the 
jet boundary of MKN501 and turbulent acceleration at the 
jet spine. However £ is usually related to the turbulent spec- 
tral index (|Biermann fc Strittmatter] (|l987l )) which may be 
different at the jet boundary and jet spine. 



5 DISCUSSION 

As the AGN jet moves through the ambient medium the 
viscosity involved will cause a shear at the jet boundary 
and hence acceleration of particles in these shear layer is 
unavoidable. If the shear gradient dT/dr is very steep or if 
the shear layer is very thin (|14fl , then shear acceleration can 
dominate over the turbulent acc elerati o n init iated by the 
instabilities at the jet boundary! Eilck (|1982j)). Turbulent 
acceleration ma y play an importan t role at the interior re- 
gions of the jet ( Virtanen fc Vainlol (120051 1 ) and can provide 
an alternative to ex plain the emissio n from the in t er kn ot 
regions of AGN jets (jMacchettd (|l996h ; IJester et al l J200ll)1 
The o bserved hard spectra at the jet spine ( Giroletti et al.l 

(2004) ) also supports this inference since turbulent accelera- 
tion c an produce a hard particle spectra (|Virtanen fc Vairiiol 

(2005) ) (also shown in section §(|3J). The electrons acceler- 
ated by the turbulence can be reaccelerated by shocks and 
can form a broken power law electron spectrum. This can 
possibly explain the break in the radio-to- x-ray spectra of 
the knots of FRI je ts JSahavanathanl l|2008h ). 

iGiroletti et al.l (|2004l ) calculated the jet viewing angle 
(8) using th e correlation between th e core power and the 
total power (|Giovannini et all l|200ll )). They estimated the 
jet viewing angle to be within 10° < 8 < 27° by compar- 
ing the observed core radio power and the expected intrin- 
sic core power derived from the correlation. However this 
estimation may vary if the core flux density variability is 
more than factor 2. Also considering the variation of the 
parameter values in the correlation with increase number 
of samples, this may not provide a strong constrain on the 
jet viewing angle. The estimate of 8 b ased on the adiaba t- 
ically expanding relativistic jet model (|Baum et al.l (|l997T 
may not be a strong constraint as it considers a simplified 
situation. Also the constrain is less severe in case of per- 
pendicular magnetic fields and observed polarisation stud- 
ies have indicated the presence of perpen d icular magnet ic 
fields at jet spine dPushkarev et aTT (120051 ); lAaronl i|l999h '). 
IStawarz fc Ostrowskil (120021 ') proposed a model similar to the 
present one, however their aim was to show the observa- 
tional implications of the two-component particle spectrum 
(power law distribution with high energy pile-up) formed at 
the boundary shear layer and the complex beaming pattern. 



6 CONCLUSION 

The observed limb brightened structure seen in the radio 
maps of MKN501 jet can be explained if we consider the 
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shear acceleration of particles at the boundary due to ve- 
locity stratification and their diffusion into the jet medium. 
This inference does not demand large viewing angle which 
is required otherwise for the explanation via differential 
Dopplcr boosting of the jet spine and boundary. We have 
shown that shear acceleration dominates over turbulent ac- 
celeration at the boundary if we consider thin shear layer or 
a sharp velocity gradient. Also for the estimated set of pa- 
rameters, shear acceleration timescale is much smaller than 
synchrotron cooling timescale allowing acceleration of elec- 
trons to be possible. The thickness of the limb brightened 
structure will be decided by the distance electrons have dif- 
fused into the jet medium before loosing its energy via syn- 
chrotron radiation. However the estimated thickness is be- 
yond the resolution of present day telescopes. Simple ana- 
lytical solution of the steady state diffusion equation consid- 
ering mono-energetic injection and particle escape, indicates 
a steep particle spectra for the electrons accelerated at the 
shear layer in comparison with turbulent acceleration. The 
radio spectral index map of MKN50f jet is also observed to 
have steep spectra at the boundary supporting the presence 
of shear acceleration. 

The author is grateful to the anonymous referee for 
his useful comments which helped in clearing many of the 
ignorances and a better understanding. The author ac- 
knowledges the useful discussions with S. Bhattacharyya, N. 
Bhatt, M.Choudhury and A.Mitra. The author is grateful to 
L. Stawarz and F. M. Rieger for enlightening information on 
various topics related to shear acceleration. 
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